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Abstract 

We theoretically study, through combining the density functional theory and an unfolding tech¬ 
nique, the electronic band structure and the charge doping effects for the deposition of potassium 
(K) on multilayer FeSe films grown on SrTiOs (001) surface. These results form a theoretical 
base line for further detailed studies of low-temperature electronic properties and their multiway 
quantum engineering of FeSe thin films. We explain the Fermi surface topology observed in exper¬ 
iment and formulate the amount of doped electrons as a function of atomic K coverage. We show 
that the atomic K deposition efficiently dopes electrons to top layer FeSe. Both checkerboard and 
pair-checkerboard antiferromagnetic (AFM) FeSe layers show electron pockets at M point and no 
Fermi pocket at F point with moderate atomic K coverage. The electron transfer from K adsorbate 
to FeSe film introduces a strong electric field, which leads to a double-Weyl cone structure at M 
point in the Brillouin zone of checkerboard-AFM FeSe. We demonstrate that with experimentally 
accessible heavy electron doping, an electron-like Fermi pocket will emerge at F point, which should 
manifest itself in modulating the high-temperature superconductivity of FeSe thin films. 
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The superconductivity in monolayer FeSe film grown on SrTiOs (001) surface has recently 
been observed [l-9|. It is attracting extensive studies on this unique system 10-21|. The 
monolayer FeSe film is the fundamental building block of crystalline iron-based supercon¬ 
ductors. The FeSe/SrTiOs system provides us with a nice platform to study the monolayer 
FeSe film thoroughly, and may shed light on superconducting mechanism of the entire iron- 
based superconductors. The superconducting transition ternperature Tc of FeSe/SrTiOa 


exceeds that of bulk FeSe by one order of magnitude 
FeSe layers on FeSe/SrTiOs 


22 


25 


24| . Whereas, further grown 


27| and FeSe layers on graphene 


and CaFo 


291 do not 


show such a high-T^ enhancement, which indicates a critical role played by the interaction 
between the monolayer FeSe and SrTiOs substrate. This interface-induced T^ enhance¬ 
ment is supposed to be closely associated with the electron doping from SrTiOs substrate 
30-34], as well as with the strength of electron-phonon interaction [l. 


35 


Angle re¬ 


solved photoemission spectroscopy (ARPES) measurements have found that there are only 
electron-like Fermi pockets that locate at M point, while the hole-like Fermi pocket at F 


point disappears |2|, [S], l32| . This peculiar Fermi surface topology indicates that the system 
is strongly electron doped. The amount of doped electrons is increased with increasing an¬ 
nealing time, and hnally reaches a saturated value. Recently, by depositing potassium (K) 
atoms onto multilayer FeSe films grown on SrTiOs surface (we henceforth name the system 
as K/n-FeSe/SrTiOs, where n is the number of FeSe layers), Miyata et al. demonstrated 
that atomic K deposition can dope electrons to FeSe Elms more effectively 40l|. In this way, 
the electron doping can be tuned continually and reversely, thus expanding the accessible 
doping range towards the heavily electron-doped region. The superconducting ^ase dia¬ 
gram shows a clear dome-like shape and the maximum T^ is as high as 48 Kelvin ^ . Along 


this line, further experimental measurements are being carried out 


41 


42|. 


To self-consistently depict the above-mentioned charge doping and predict new adsorbate- 
tuned electronic properties superposed onto the inherently profound physics of the 
FeSe/SrTiOs system, apparently, a systematic and enlightening theoretical study of K de¬ 
position affection to the electronic structures of FeSe films is essential but still keeps to 
be initialized. In this paper, by using the density functional theory (DFT) based hrst- 
principles calculations and checkerboard (CB) as well as pair-checkerboard (PCB) antifer¬ 
romagnetic (AFM) models, the energy band structure and the charge transfer effects of the 
K/n-FeSe/SrTiOa system is studied. The Fermi surface topology observed in experiment is 
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FIG. 1: (Color online) The side (a) and top (c) view of K/2-FeSe/STO atomic structure. Panel 
(b) shows the two-dimensional hrst Brillouin zone of K/n-FeSe/STO. 

satisfactorily explained. We find that only the top layer FeSe is effectively charged. The 
amonnt of doped electrons linearly depends on atomic K coverage (denoted by px)- The 
induced electric held strongly distorts the CB FeSe energy bands, leading to a double-Weyl 
cone structure at M point. Remarkably, our results also show that an electron-like Fermi 
pocket will rapidly develop at F point if the amount of doped electrons further increases to 
be larger than 0.2 e per FeSe formula unit in CB FeSe. For PCB FeSe, the critical value is 
0.17 e. 

The DFT calculations in this work were mainly performed by using Vienna Ab-initio 


Software Package (VASP) |43|. T 


le atom core electrons were described by projector aug¬ 


mented wave (PAW) method 4^, l45| . Perdew-Burke-Ernzerhof 46[ (PBE) functional was 
used to treat the electronic exchange-correlation. The energy cutoff for the plane-wave basis 
was set at 400 eV. The hrst Brillouin zone was sampled in the fc-space with Monkhorst-Pack 
scheme and grid sizes are 13x13, 7x7, 5x5 and 8x4 for 1x1, 2x2, 3x3 and 2x4 supercells, 
respectively. We have checked that the total energy is converged for the cutoh energy and 
A;-point sampling. The atomic structure was relaxed until the force on each atom is smaller 
than 0.01 eV/A. The K/n-FeSe/SrTiOs can be divided into three parts (K atoms, FeSe layers 
and SrTiOs substrate). Considering that weak attractions exist in this composite system. 


we have added the van der Waals correction to DFT calculations. The slab models were 
used and more than 12 A vacuum space was added to avoid artihcial interaction. The in¬ 


plane lattice parameter was set to 3.901 A 3^. We used 6-layer SrTiOa with the Ti02 layer 
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FIG. 2: (Color online) The energy band structures of K/2-FeSe/STO (a) and K/FeSe (b) isolated 
systems. In panel (a), the energy bands projected to K atom and top layer FeSe are signed by blue 
dots. The Fermi level is set to be zero. Panel (c) shows the charge transfer from K atom to the 
top few FeSe layers in a K/n-FeSe system. The inner hgure shows the K/n-FeSe atomic structure. 


on the top to mimic the SrTiOs substrate. Besides the VASP code, we also performed DFT 
calculations in Quantum Espresso (QE) 48(] code with norm-conserving pseudopotentials 


49| to get accurate unfolded electronic energy bands by using a home-built code. 


The relaxed atomic structure of the K/ 2 -FeSe/SrTi 03 is shown in Figs. [I](a,c). The px is 
1, which is one K atom per FeSe primary cell. The thickness of the down layer FeSe, which 
is dehned by the distance between the top Se atom in FeSe and top layer Ti02, is 5.62 A. 
The thickness of the top layer FeSe is 5.47 A, which is very close to the c lattice parameter 
of bulk FeSe 2^. The K atoms adsorb on top layer FeSe, right on the top of bottom Se 
atoms. The distances between the K atom and top Se atoms and Fe plane are determined 
to be 3.24 A and 3.14 A, respectively. 

The energy bands of K/ 2 -FeSe/SrTi 03 with CB magnetic orientation are shown in Fig. 
[2](a). All the energy bands are plotted in green color, and the energy bands that are mainly 
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contributed by top layer FeSe and K atom are signed by blue dots. To compare with these 
energy bands, we also calculated the energy bands of K/FeSe system, which consists only of 
the K atoms and top layer FeSe. The resultant band structure is shown in Fig. [2]^b). It is 
quite similar to the blue dots in Fig. Wisi), while the energy bands in Fig. Ula) without blue 
dots are quite similar to those in Fig. 2(a) in Ref. 3^. This means that the interaction 
between FeSe layers is quite weak, and the atomic K adsorption only affects the top layer 
FeSe. We then further checked the energy bands of K/S-FeSe/SrTiOa system with CB 
magnetic orientation and K/ 2 -FeSe/SrTi 03 with PCB magnetic orientation. The results 
are the same by the observed fact that the adsorbed K atoms only influence the top layer 
FeSe and have negligible effects on the other FeSe layers. The doping charge distribution in 
a multilayer FeSe system can be calculated by using a K/n-FeSe slab model. The calculation 
results of K/5-FeSe model with different px and magnetic orientations are shown in Fig. 
|2](c). All the calculation results show that the top layer FeSe has much more doping electrons 
than the second layer (more than 31 times larger). The doping electrons in the non-top layer 
FeSe is negligible. This phenomenon agrees with the most recent experimental results 50|. 
In other words, the K/FeSe slab model can approximately describe the K adsorption on 
n-FeSe/SrTiOa (n> 2) systems correctly. Thus, in this work, we studied the detailed K 
adsorption effects by using K/FeSe slab model with different values of px- 

We then studied the energy bands of CB FeSe with different px- The energy bands 
can be calculated by using a primary cell, if the is 0 or 1. The corresponding results 
are shown in Figs. I3i:a,e). For the cases of fractional px, supercells have to be used to 
calculate their electronic structures. However, the supercell energy bands are quite different 
from primary cell energy bands. In order to compare energy bands between the primary 
cell and the supercell, we resort to unfolding the supercell energy bands by using a home- 
built code In PAW calculations, the norm of PAW wavefunction is not restricted to 


be unity. Therefore, the energy band unfolding weight within PAW is not accurate 511. 
Thus, we performed the QE calculations with norm-conserving pseudopotentials to unfold 
the supercell energy bands accurately. The plane wave cutoff energy was 1490 eV. The 
/c-space samplings were set the same with VASP calculations. We tested that the FeSe and 
K/FeSe primary cell energy bands calculated by QE agree well with VASP results. After 


the QE supercell calculations, the WannierQO code 


52l | was used to calculate the maximum 


localized Wannier functions (MLWFs). We added a few Fortran lines to WannierQO to write 
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out the MLWFs in a format that is easier for further calculations. After that, our Quantum 


Unfolding code 


51 


53| reads these MLWFs and unfold the electronic energy bands to the 


hrst Brillouin zone of FeSe primary cell. 

We used 2x2 CB FeSe supercells with 1, 2 and 3 K atoms adsorbed on one side of FeSe 
monolayer to mimic the K/n-FeSe/SrTiOa systems. The corresponding px are 0.25, 0.5 and 
0.75 respectively. The unfolded energy bands are shown in Figs. |3]^b-d). For comparison, 
the energy bands for the cases of pk=0 and pk=^ (by usage of 1x1 primary cells) are also 
displayed in Fig. |3] [panels (a) and (e)]. Initially [panel (a)], the Fermi surface contains two 
parts. They are a small electron-like Fermi pocket at M point and a hole-like Fermi pocket 
at F point. With increasing p^ step by step, then, the energy bands continually go down 
with respect to the Fermi level. As a consequence, the hole-like Fermi pocket disappears, 
while the electron pocket enlarges its size dramatically. Remarkably, these Fermi surface 
transformations shown in Figs. Ha-c) agree well with the most recent ARPES experiments 
40(1, in which the K depositions with the electron doping up to 0.15 e per FeSe formula 
unit were explored. Prominently, as we increase px, a conduction band (indicated by a 
green arrow) at P point lowers its energy much faster than other energy bands. It closes 
to the Fermi level at px = 0.5, and crosses the Fermi level at px = 0.75. The crossing 
introduces an electron-like Fermi pocket at F point. This agrees with the previous energy 
band calculation in uniformly electron doped FeSe monolayer [^. Our unfolded energy 
bands show that this phenomenon is robust against the doping charge distribution, and may 
be realized in K/n-FeSe/SrTiOa system. It is a Lifshitz transition since the Fermi surface 
topology is altered. The effect of this Lifshitz transition to FeSe monolayer superconducting 
is still not clear, therefore need more theoretical and experimental study. 

On the other hand, comparing to the case of uniformly electron doped FeSe, the electron 
doping raised by atomic K adsorption brings about a new feature. That is, the energy bands 
at M point become non-symmetric, as shown by the red arrows in Fig. 3. They move to the 
right side (M-Y direction). The bias direction depends on the electron spin. We have checked 
that the spin-down energy bands bias to the left side (M-X direction). This phenomenon 
also exists in l-FeSe/SrTiOs system j^. Our DFT and tight-binding calculations on FeSe 
monolayer in electric held revealed this peculiar feature, and concluded that the energy bands 
biased at M point is caused by a strong electric held along the perpendicular direction of 
FeSe him. It is a self-established electric held as a consequence of charge transfer. The 
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FIG. 3; (Color online) Spin-up unfolded electronic energy bands of checkerboard-AFM FeSe mono- 
layer with pk = 0.25 (b), 0.5 (c) and 0.75 (d). The spin-up electronic energy bands of FeSe primary 
cell with (e) and without (a) K atom adsorption are also shown here. Panel (f) shows the warp 
gap at M point and electric potential difference AF as functions of pK- Panels (g) and (h) show 
the double-Weyl cones around M point for up and down spin electrons respectively. 

strong electric field leads to an electric potential difference between top and bottom Se 
atoms (AF=VJ;op-Se—f4ot-Se)- The top (bottom) Se atoms are at the ±X (±F) directions 
of a Fe atom with spin-up magnetic moment, and are at ±Y (±A) directions of a Fe atom 
with spin-down magnetic moment. Thus, AV will lead to a difference between X and Y 
directions. We studied the energy bands around the red arrow and hnd that they compose 
a double-Weyl cone structure. Different from Dirac cone, the Weyl cone does not have spin 
degeneracy. The sketches of double-Weyl cones for spin-up and spin-down electrons are 
shown in Figs. |3](g,h). The two Weyl points for up spin are in M-Y direction, while the two 
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Weyl points for down spin are in M-X direction. We name the energy gap at M point as 
a warp gap, as shown in Fig. |3t^d). It characterizes the warping strength of energy bands. 
The calculated warp gap, as shown by the red line in Fig. El^f), is found to linearly increase 
with the K coverage at the range of 0 < px < 0.5. Its slope decreases when px > 0.5. We 
ooint out that such kinds of warp gaps have been observed in the l-FeSe/SrTiOs system 
TJ. The htted yellow line in Ref. 7| indicates the upper part of the double-Weyl cones in 
M-F direction in Fig. 1(b) of Ref. 7|], wherein the bottom of the yellow line clearly departs 
from the lower energy bands. The estimated warp gap is about 50 meV. The wa rp g ap can 
also be estimated from other ARPES experiments in l-FeSe/SrTiOs system 25|, l32|. From 
Figs. 2(b,c) in Ref. |32| and Fig. 2(d) in Ref. |25|, we estimated that the warp gap is 
about 50 meV too. Based on our calculations, to observe the double-Weyl cone energy band 
structure, the ARPES spectra should be carried out along the M-X and M-Y directions. The 
ARPES spectrum would show four Weyl cones around M point both along M-X and M-Y 
directions, if the measurement does not distinguish the up and down electron spins. From 
previous studies [S^, we know that the warp gap relates to AY directly. Thus, we calculated 
the value of AV too. The electronic potential at a Se atom was obtained by averaging the 
electronic potential at atom core area. The calculated AY is shown by the blue line in Fig. 
[3](f). It linearly depends on px in the whole K coverage range we considered. 

The unfolded PCB K/FeSe energy bands are shown in Figs. |l](a-d). For px = ^ case, 
the unfolded energy bands in Fig. |4l(a) show that there are double-Dirac cones around M 
point. It agrees with ARPES experiment results. Unlike the double-Weyl cones in CB 
FeSe, the double-Dirac cones in PCB FeSe exist without doping electrons. The connecting 
line between the two Dirac points is in F-M direction. They are the same for spin-up 
and spin-down electrons, as shown by the sketches in Figs. Ill(e,f). The unfolding energy 
bands for px = ^ show that there is no Fermi surface exists at F point. The highest 
valence energy bands locates at about 100 meV below Fermi energy. This conflicts with 
ARPES experiments, which show the existence of hole-like Fermi pocket at F piont. After 
K atom deposition (px=0.25), the whole energy bands move downwards with respect to 
Fermi energy. The Dirac points locates about 200 meV below Fermi energy. And the Fermi 
surfaces around M point combines to form a circle, as shown in Fig. 111(b). The PCB FeSe 
Fermi surface topology at px=0A5 is the same as that of CB FeSe. And it is consistent with 
ARPES experiment results. After we increase K atom coverage to px=0A, an electron-like 
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FIG. 4: (Color online) Spin-up unfolded electronic energy bands of PCB FeSe monolayer with 
Pk = 0 (a), 0.25 (b), 0.5 (c) and 0.75 (d). Panel (e) and (f) show the double-Dirac cones around 
M point for spin-up and spin-down electrons, respectively. 


Fermi pocket exists at F point. It originates from the energy band bottom signed by a green 
arrow in Figs. |ll(a-c). The energy band bottom crosses the Fermi level at 0.25< px <0.5. 
The Pk at the crossing point is smaller than that in CB FeSe, which is 0.5< px <0.75. As 
we further increase K atom coverage to pi^=0.75, the energy band bottom at F point lowers 
its energy, which increases the size of electron-like Fermi pocket. 

From the energy band variations, we know that atomic K adsorption efficiently dopes 
electrons to the top FeSe layer. At this part, we examined the amount of doped electrons 
by performing Bader analysis 5^. In PCB FeSe, we used the \/2 x 2\/2 supercell models 
with 0~4 K atoms. The calculation results are shown in Fig. [5] by green blocks. In CB 
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FeSe, besides primary cell model for FeSe and K/FeSe, we also used three kinds of supercell 
models. They are 2x2 supercell with 1~3 K atoms, \/E x \/E supercell with 1~4 K atoms 
(Pk = 0.2, 0.4, 0.6, 0.8), and 3x3 supercell with 1 K atom {px = ^)- The Bader analysis 
results are shown in Fig. [5] by red down triangles. Besides the PCB and CB phases of 
FeSe monolayer, we also studied the K adsorption on paramagnetic phase (PM) of FeSe 
monolayer. The calculated electron transfer results are shown in Fig. [5] by up triangles. 
The electron transfers for PCB, CB and PM FeSe monolayers are the same for px < 0.5, 
and have a slight difference for px > 0.6. The maximum deviation is 0.016 e, locating at 
Px = 0.75. Thus, we conclude that the magnetic state of FeSe has negligible effect on the 
electron transfer from the adsorbed K atoms. The electron doping increases with increasing 
Px- It is linear when px < 0.5. We htted calculated data points (in unit of e) at this range 
by using Uf, = c x px- The rig here is the amount of doped electrons per FeSe formula unit. 
The htted slope c is 0.42. The value of rie deviates from the Ue = c x px line for px > 0.5. 
It is lower than the linear line value. The knee point locates between px = 0.5 and 0.6. It 
is near the place where Lifshitz transition may happen in PCB and CB FeSe layer. 

The reported maximum value of electron transfer in experiment j^] is 0.15 e. From 
the linear relation Ug = c x px, we estimate that the px is about 0.36. It is still in the 
linear range. The corresponding CB FeSe energy band structure is in the interval between 
those shown in Figs. [3](b) and [3](c), where the electron pocket at F point still does not 
emerge. The PCB FeSe energy bands for px = 0.36 is in the interval between those shown 
in Figs. 111(b) and 111(c). It is the same interval where Lifshiz transition occurs. However, our 
linear extrapolation show that Lifshiz transition in PCB FeSe occurs at px = 0.37, which 
is larger than 0.36. Thus, electron pocket at F point does not emerge. Both the calculation 
results at PCB and CB FeSe are consistent with the ARPES experiment 40|, which did not 
hnd any Fermi surface around F point. From the energy band structures in Figs. l3](a-e) 
and Figs. Ha-c), and the linear relation rig = c x px, we estimate that the electron-like 
Fermi pocket at F point would appear for Ue > 0.2 e in CB FeSe and Ug > 0.17 e in PCB 
FeSe. To achieve this, the simplest way is to deposit more K atoms on top layer FeSe. To 
avoid the potential K-clustering problem sometimes occurred at large coverages, here, we 
propose an alternative experiment that can dope more electrons to FeSe monolayer. It is 
K deposition on monolayer FeSe grown on SrTiOs (001) surface. After sufficient annealing, 
SrTiOs substrate dope some amount of electrons to FeSe monolayer. Then the K deposition 
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FIG. 5: (Color online) The charge doping (ng) as a function of K atom coverage {pk) in paramag¬ 
netic (green blocks), checkerboard-AFM (red down triangles), and pair-checkerboard-AFM (blue 
up triangles) FeSe. The dotted black line is fitted from small px data. 


on FeSe monolayer will fnrther increase the electron doping. The maximnm Uf. is probably 
larger than 0.2 e. Interestingly, a previons experiment has alrea^ found an electron energy 
band centered at F point, existing 75 meV above Fermi energy |9|. It is very close to Fermi 
energy. The energy bands in this experiment would lie in the interval of Figs. |3](b) and 
He), and is similar to Fig. Hb)- Thus, we estimate that once we dope more electrons to 
the system by K deposition, the F centered energy band would cross the Fermi level and 
introduce an additional electron-like Fermi pocket. 

In summary, we have theoretically studied the atomic K deposition on n-FeSe/SrTiOs 
system. By using CB and PCB FeSe models, we explained the Fermi surface topology 
observed in experiment, i.e., the existence of electron-like pockets at M point, as well as the 
disappearance of the hole-like pockets and the enlargement of the electron-like pockets after 
K deposition. Our results show that the K adsorbate only affects the top layer FeSe. The 
charge transfer linearly depends on the K coverage in the range 0 < < 0.5 and leads to a 

strong electric field along the perpendicular direction of FeSe film. This strong electric field 
warps the CB FeSe energy band structure at M point and causes the occurrence of a double- 
Weyl cone structure. These Weyl points locate at M-Y and M-X high-symmetry lines for 
spin-up and spin-down electrons, respectively. Also, our energy band unfolding calculation 
shows that PCB FeSe contain double-Dirac cones at Femi energy around M point for px = 
0. They are along F-M direction. We have also shown that an electron-like Fermi pocket will 
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emerge at T point for Ue > 0.2 e in CB FeSe and Ue > 0.17 e in PCB FeSe, which is expected 
to be observed in the undergoing atomic K deposition experiments. Our obtained double- 
Weyl cones at M point and electron-like pocket at F point are typical features of CB FeSe. 
The spin-resolved experiments in the future will hgure out whether CB or PCB magnetic 
orientation is the basement of the subsequent profound physics for the FeSe monolayer and 
thus increase our knowledge about this unique low-dimensional material. 

This work was supported by Natural Science Foundation of China (NSFC) under Grant 
Nos. 11474030 and 91321103, and by the National Basic Research Program of China (973 
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